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Abstract 

Single-crystalline CeRh3Si2 was investigated by means of x-ray diffraction, magnetic suscep- 
tibility, magnetization, electrical resistivity, and specific heat measurements carried out in wide 
temperature and magnetic field ranges. Moreover, the electronic structure of the compound was 
studied at room temperature by cerium core- level x-ray photoemission spectroscopy (XPS). The 
physical properties were analyzed in terms of crystalline electric field and compared with results 
of ab-initio band structure calculations performed within the density functional theory approach. 
The compound was found to crystallize in the orthorhombic unit cell of the ErRh3Si2 type (space 
group Imma - No. 74, Pearson symbol: oI24) with the lattice parameters: a = 7.1330(14) A, 
b = 9.7340(19) A, and c = 5.6040(11) A. Analysis of the magnetic and XPS data revealed the 
presence of well localized magnetic moments of trivalent cerium ions. All physical properties were 
found to be highly anisotropic over the whole temperature range studied, and influenced by ex- 
ceptionally strong crystalline electric field with the overall splitting of the 4/^ ground multiplet 
exceeding 5700 K. Antiferromagnetic order of the cerium magnetic moments at Tn = 4.70(1) K 
and their subsequent spin rearrangement at Tt = 4.48(1) K manifest themselves as distinct anoma- 
lies in the temperature characteristics of all investigated physical properties and exhibit complex 
evolution in an external magnetic field. A tentative magnetic B—T phase diagram, constructed for 
B parallel to the 6-axis being the easy magnetization direction, shows very complex magnetic be- 
havior of CeRh3Si2, similar to that recently reported for an isostructural compound CeIr3Si2. The 
electronic band structure calculations corroborated the antiferromagnetic ordering of the cerium 
magnetic moments and well reproduced the experimental XPS valence band spectrum. 

PACS numbers: 75.10.Dg (crystal-field theory and spin Hamiltonians) . 75.30.Gw (magnetic anisotropy), 
75.30.Kz (magnetic transitions), 75.50.Ee (antiferromagnetics), 71.20.Lp (electronic structure of bulk mate- 
rials) 
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I. INTRODUCTION 



The phase diagram Ce-Rh-Si comprises twenty five ternary phases with well defined crys- 
tal structures l|. However, physical properties have been reported for only a few of them, 
namely: CeRhSi2 and Ce2Rh3Si5 (both described in the literature as intermediate valence 



systems 



and CeRhSis (antiferromagnetic heavy-fermion systems (6|j7|). The two latter phases are 



known as pressure-induced superconductors [SHllL While in CeRh2Si2 the superconduc 



Ce3Rh3Si2 (exhibiting complex long-range magnetic ordering jsl), CeRh2Si2 



the ordered state 
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tivity seems to replace the antiferromagnetism js Q], in CeRhSia it develops deeply inside 



Most interesting, the superconductivity in non-centrosymmetric 



CeRhSis involves spin-fiuctuations-assisted triplet pairing In the context of recent find- 
ings for the Ce-Rh-Si ternaries, comprehensive physical characterization of other (hitherto 
unknown or hardly studied) phases from the system seems quite desirable. To the best of 
our knowledge existence of CeRh3Si2 has for the first time been reported by Morozkin et al. 



12|, yet without any crystallographic or physical data. Then, CeRh3Si2 has been found as 



an impurity phase in polycrystalline CeRh2Si2 [13[. In the latter study it has been estab- 
lished that the compound crystallizes with an orthorhombic crystal structure and exhibits 
an antiferromagnetic phase transition at Tn = 5 K [l^. 

Here we report the results of our investigation of the crystal structure and the basic 
physical properties of CeRh3Si2 performed on single-crystalline specimens in wide temper- 
ature and magnetic field ranges. The experimental characterization is accompanied by the 
results of ab-initio calculations of the electronic band structure. Short accounts on the work 
presented in this paper were given in conference communications 



II. EXPERIMENTAL AND COMPUTATIONAL DETAILS 

A single crystal of CeRh3Si2 was grown by the Czochralski pulling method employing a 
tetra-arc furnace under protective ultra-pure (4.8N, BOC Gases) argon atmosphere, which 
was additionally titanium-gettered during the whole growing process. The starting poly- 
crystalline melt was prepared from pieces of Ce (3N, Ames Laboratory), Rh ingot (3N, 
Chempur) and Si chips (6N, Chempur). The pulling rate was 10 mm/h, and the copper 
heart rotation speed was 3 rpm. The final ingot of CeRh3Si2 was about 4 mm in diameter 
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and 40 mm in length. From this single-crystaUine rod a fragment of about 20 mm in length 
was cut for experimental studies. The crystal was wrapped in Ta-foil, sealed in an evacu- 
ated silica tube, and annealed at 900°C for 2 weeks. The quality of the product was verified 
by means of x-ray powder diffraction measurements performed on a pulverized piece of the 
crystal (Stoe diffractometer with Cu Ka radiation) and microprobe analysis (Phillips 515 
scanning electron microscope equipped with an EDAX PV 9800 spectrometer). The entire 
X-ray diffractogram was indexed with an orthorhombic unit cell, and the sample was found 
to be single phase. Additionally, a polycrystalline sample of an isostructural LaRh3Si2 com- 
pound, being a non-magnetic counterpart of CeRh3Si2, was prepared by conventional arc 
melting stoichiometric amounts of the constituents. X-ray powder diffraction experiments 
proved that the La-based phase is indeed isostructural to CeRh3Si2. 

The crystal structure of CeRh3Si2 was also examined on a small single-crystalline frag- 
ment using a four-circle diffractometer equipped with a CCD camera (Kuma Diffraction 
KM-4 with graphite-monochromatized Mo Ka radiation). Crystal structure refinement was 
performed using full-matrix least-squares on as a refinement method and applying the 
SHELXL-97 program 16|]. The magnetic properties were studied at temperatures ranging 
from 1.7 K up to 300 K in applied magnetic fields up to 5 T using a Quantum Design 
MPMS-5 superconducting quantum interference device (SQUID) magnetometer. The elec- 
trical resistivity was measured in the temperature range 1.5 - 300 K employing standard 
DC four-point technique. Electrical contacts to bar-shaped specimens were made by spot 
welding. The heat capacity was measured in the temperature interval 350 mK - 300 K in 
applied magnetic fields up to 9 T using a Quantum Design PPMS platform. For these ex- 
periments the thermal relaxation method 17| was applied with 2% heat pulses and Apiezon 
N vacuum grease as a sample mounting medium. The x-ray photoelectron spectroscopy 
(XPS) measurements were carried out at room temperature with monochromatized Al Ka 
radiation using a PHI-5700 ESCA spectrometer. The sample of CeRh3Si2 was scraped with 
a diamond file under high vacuum immediately before recording a spectrum. 

The electronic structure calculations were performed with a density functional theory 



(DFT) [l8j using the full-potential linearized augmented plane wave (FP LAPW) method 
implemented in the latest version (Wien2k) of the original WIEN code 19|. In the Wien2k 
code calculations the scalar relativistic approach was implemented with the spin-orbit in- 
teractions taken into account by using the second variational method 2^. The exchange 
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correlation potential in the local (spin) density a ppr oximation (LSDA) was assumed in the 



form proposed by Perdew, Burke, and Ernzerhof [2]| using the generalized gradient approx- 
imation. Furthermore, to improve the description of the strongly correlated 5f electrons, the 
on-site Coulomb energy U correction and exchange J parameters were introduced within 
the LSDA+U approach [22|j23|. In our case we used Ucs = U — J, setting J = 0. Another 
approach applied by us to reduce the discrepancy between the DFT magnetic moments and 
experimental ones was taking into account the so-called orbital polarization (OP) term, as 



proposed by Brooks 



and Eriksson et al. 



25 



26j. The number of fc-points was 2000 in 



the Brillouin zone (BZ), which corresponds to 240 points in an irreducible wedge of the BZ 
for all methods of calculations applied in this paper. For the BZ integration, a tetrahedron 



method was used 



271 ]. The self-consistency criterion was equal 10 ^ Ry for the total energy. 



The calculations were performed for lattice constants and atomic positions in the unit cell 
described below in Sec. IIII A[ The theoretical x-ray photoemission spectra were obtained 
from the calculated densities of electronic states (DOS) convoluted by a Gaussian with a 
half- width S equal to 0.3 eV and scaled using the proper photoelectronic cross sections for 
partial states [28 1. 



III. RESULTS AND DISCUSSION 



A. Crystal structure 

The crystal structure refinement (see Tab. [T]) performed for single- crystalline CeRh3Si2 
confirmed that the compound crystallizes in the orthorhombic symmetry with the space 
group Imma (No. 74, Pearson symbol oI24) and that it is isostructural with ErRh3Si2 29 |. 
Refined lattice parameters are: a = 7.1330(14) A, b = 9.7340(19) A, c = 5.6040(11) A, and 
the atomic coordinates are listed in Tab. HIl Figure [H presents the structure of CeRh3Si2 
viewed along the a-axis. 



The ErRh3Si2-type crystal structure was described in detail in Ref. 
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It represents a 

deformed superstructure of the hexagonal CaCus-type structure with a = 2ch, b = V^a^, 
and c = a/j, where ah and Ch are the lattice parameters of the parent hexagonal unit cell. 
The characteristic features of this atom arrangement are a honeycomb- like network of the 
Rh and Si atoms within the be plane and zigzag chains of the Ce atoms along the a-axis (see 
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FIG. 1: (Color online) Crystal structure of CeRh3Si2. 



Fig. [I]). 

B. Magnetic properties 

Figures [2](a,b) present the magnetic susceptibility x of CeRh3Si2 measured as a function 
of temperature T in a magnetic field B of 0.5 T applied parallel to the main crystallographic 
directions. Sharp susceptibility maxima at about 4.5 K, best visible on the Xb{T) curve 
measured along the 6-axis [see Fig. [2](a)], manifest the onset of an antiferromagnetically or- 
dered state. Below 4.5 K the susceptibility Xb does not decrease significantly, as expected for 
simple antiferromagnets, but saturates at low temperatures, hence suggesting more complex 
magnetic behavior. The anomaly at 4.5 K is only slightly weaker for the c-axis, and the 
behavior of Xc in the ordered region is similar to that of Xb- On the contrary, the a- axis 
susceptibility is featureless and very small, thus indicating that the magnetic moments are 
confined to the be plane of the unit cell. A distinct difference of XaiT) as compared to the 
other two components is observed also in the paramagnetic state [see Fig. E^b)]. At room 
temperature, the ratio Xb/Xa is about 4, and it increases with decreasing temperature down 
to 200 K, where it is larger than 6. The Xa^i'^) variation exhibits strongly curvilinear shape 
in the whole T-range studied. In contrast, the x^^(T) and x~^(T) curves are almost identi- 
cal and exhibit nearly linear behavior characteristic of systems with well localized magnetic 
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FIG. 2: (Color online) (a) Low-temperature magnetic susceptibility of single-crystalline CeRh3Si2 

measured in external magnetic field B applied parallel to the principal crystallographic axes. The 

arrows mark the Neel temperature; the solid lines serve as a guide for the eye. (b) Temperature 

variations of the inverse magnetic susceptibility of the compound; x^vg represents the inverse 

average magnetic susceptibility. The solid line is a fit of the Curie- Weiss law to the experimental 

data; the dashed curves represent the fits described in Sec. IIIIEI (c) Magnetization as a function 

of applied magnetic field measured for single-crystalline CeRh3Si2 in zero- field (bright symbols) 
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and field cooling (dark symbols) regimes, along the principal crystallographic axes. The solid lines 
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FIG. 3: (Color online) Temperature dependencies of the magnetic susceptibility of single-crystalline 
CeRh3Si2 measured in different applied magnetic fields B parallel to the 6- axis. For the sake of 
clarity the curves are shifted up. 

moments. The overall behavior of the magnetic susceptibility of CeRh3Si2 implies huge 
magnetocrystalline anisotropy due to strong crystal field interactions (cf. Sec. IIIIE|) . 

As can be seen in Fig. [2](b), the inverse average "polycrystalline" magnetic susceptibility 
X^g [defined as Xavg = |(Xa + Xb + Xc)] follows above about 150 K the Curie- Weiss law 
with the effective magnetic moment fics = 2.47 /ie and the paramagnetic Curie temperature 
^cw = —52 K [see the straight solid line in Fig. [21^b)]. The experimental value of fies is close 
to that calculated for a free Ce^"*" ion {g>y j{j + 1) = 2.54), and thus indicates the presence 
of well localized magnetic moments. The negative value of 6'cw hints at antiferromagnetic 
character of the magnetic exchange interactions, in line with the AFM ordering observed 
below 4.5 K. 

Figure[2](c) shows the magnetization of CeRh3Si2 measured at 1.72 K as a function of the 
magnetic field strength. Here the complex magnetic behavior of the compound manifests 
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FIG. 4: (Color online) Isothermal magnetization of single-crystalline CeRli3Si2 measured with the 
magnetic field applied along the 6-axis at several different temperatures from the ordered region. 
The curves are shifted up for the sake of clarity. 

itself as cascades of metamagnetic-like transitions in M{B) measured along the b- and c- 
axes. For the former direction of the magnetic field one observes some saturation above 
about 1.5 T at a value /ig = 12.8 emu/g corresponding to about 1.16 /ib- In the case of 
B II c-axis, the magnetization increases linearly with increasing B reaching in 5 T a value of 
about 0.9 /iB- In turn, for B \\ a-axis the magnetization is a linear function of the apphed 
field in the entire range studied, and the value of M measured at 5 T is nearly 20 times 
smaller than that observed for B \\ 6-axis. The overall behavior of M{B) is reminiscent 



of that reported for CeSb |30lj3l| . TbNi2Si2 [3^, and TbNi2Ge2 [33||, which are commonly 
called "devil's staircase" systems (for short review and discussion see Sec. llVp . 

In order to shed more light on the nontypical temperature variation of the magnetic 
susceptibility below 4.5 K [Fig. EJj^a)] we performed a series of additional measurements of 
Xh(T), with increased density of the data points. As can be inferred from Fig. [3l in the 
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TABLE I: Crystallographic and structure refinement data for single-crystalline CeRh3Si2. 



Empirical formula 
Formula weight 
Temperature 
Wavelength 

Crystal system, space group 
Unit cell dimensions 

Unit cell volume 
Z 

Calculated density 
Absorption coefficient 
F(OOO) 

Q range for data collection 
Limiting indices 

Reflections collected / unique 

Completeness to ^ = 45.28 

Refinement method 

Data / restraints / parameters 

Goodness-of-fit on 

Final R indices [/ > 2cj(/)] 

R indices (all data) 

Extinction coefficient 

Largest diff. peak and hole 



Ceo.50 Rhi.50 Si 
252.51 g/mol 
293(2) K 
0.71073 A 

orthorhombic, Imma 

a = 7.1330(14) A, a = 90° 

b = 9.7340(19) A, (3 = 90° 

c = 5.6040(11) A, 7 = 90° 

389.10(13) A3 

4 

8.621 Mg/m^ 
24.354 mm-i 
884 

4.19° ^45.28° 
-13 < /i < 9 
-12 < A; < 19 
-10 < / < 10 

4366 / 840 [R{mt) = 0.0465] 
93.8 % 

Full-matrix least-squares on 

840 / / 22 

1.010 

Ri = 0.0287, wR2 = 0.0509 
Ri = 0.0432, wR2 = 0.0524 
0.0356(8) 

2.707 eA-3 and -3.148 eA'^ 
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TABLE II: Atomic coordinates and equivalent isotropic thermal displacement parameters for 


CeRhgSis. 


Ucq was defined as 


one third of the trace of the orthogonalized Uij 


tensor. 


Atom 


X [xlO^] 




y [xlO^] 


z [xlO^] 


Ucq [xlO^A^] 


Ce 







7500 


2160(1) 


8(1) 


Rh(l) 


2172(1) 




5000 





6(1) 


Rh(2) 


-2500 




7500 


-2500 


8(1) 


Si 







9215(2) 


6996(2) 


7(1) 


TABLE III: Anisotropic thermal displacement parameters ( x 


10^ A^) for CeRh3Si2. The anisotropic 


displacement factor exponent takes the form: — 27r^[/i^a^?7ii 


+ ... + 2hkabUi2] 




Atom 


Uu 


U22 


U33 


U23 


Ul3 Ui2 


Ce 


8(1) 


7(1) 


8(1) 








Rh(l) 


7(1) 


6(1) 


5(1) 


1(1) 





Rh(2) 


13(1) 


4(1) 


7(1) 





-2(1) 


Si 


6(1) 


9(1) 


7(1) 


-1(1) 






lowest magnetic field studied (i.e. B = 0.005 T) the magnetic order in CeRh3Si2 sets in 
already at Tn = 4.7 K. The Neel temperature decreases with increasing field, as expected 
for antiferromagnets. In a field of 0.5 T, Tn is reduced to about 4.5 K, marked by an 
arrow in Fig. [2t^a). Moreover, the high resolution measurements revealed the presence of 
some spin rearrangement, which takes place just below Tn and therefore is hardly visible in 
Fig. [2t^a). In 0.005 T, this spontaneous change of the magnetic structure manifests itself as 
a tiny anomaly in the x{T) curve at Tt = 4.5 K, which develops with increasing field into a 
clear kink moving down on the temperature scale. The latter observation suggests that the 
antiferromagnetic character of the magnetic order is retained below Tt. The plateau visible 
in 0.5 T below 4 K [Fig. EJ^a)] apparently results from superposition of low-temperature 
slopes of the anomalies at Tn and Tt, and another, field-induced phase transition. It is well 
visible in Fig. |3] in the field range 0.7 T ^ B ^ 0.9 T as a broad ferromagnetic-like hump 
dominating the susceptibility below Tt. Unfortunately, complex shape of the x{T) curves 
does not allow us to determine the exact temperature at which the latter contribution sets 
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FIG. 5: (Color online) Electrical resistivity of single-crystalline CeRh3Si2 measured as a function of 
temperature with electrical current i flowing along the principal crystallographic directions; solid 
lines serve as a guide for the eye. The inset: low temperature region; the arrows mark magnetic 
transition temperatures. 



Figure m presents a series of dense magnetization curves measured along the easy 6-axis. 
For the sake of straightforward comparison with the afore-discussed susceptibihty data, these 
experiments were focused on magnetic fields lower than 2 T and temperatures close to TN,f 
As seen, the "devilish" field dependence of the magnetization, observed at 1.7 K [Fig. [2](c)], 
is hardly visible at temperatures higher than 4.4 K. However, upon decreasing T a sequence 
of distinct steps in the M{B) curve emerges from the smooth curve. Three of these steps, 
located at 3.8 K at about 0.15, 0.32 and 0.53 T, dominate the magnetization, yet several less 
visible anomalies can also be found on the M{B) curve. The experimental data displayed 
on Figs. [3] and m were used for constructing a tentative phase diagram described in Sec. IIVI 

C. Electrical resistivity 

Figure |5] presents the electrical resistivity of single- crystalline CeRh3Si2 as a function of 
temperature, measured with the current i flowing along the main crystallographic directions. 
Clearly, the compound exhibits good metallic conductivity in the entire temperature range 
studied. The antiferromagnetic phase transition manifests itself as a drop in p{T) at about 



m. 
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4.7 K, being equal to the Neel temperature deduced from the magnetic measurements per- 
formed in the weakest magnetic fields (cf. Fig. |3]). Closer look at the low-temperature range 
[see the inset to Fig. unambiguously shows the presence of another kink on the p(T) curve 
that occurs at about 4.5 K. The latter finding is in-line with the scenario of two subsequent 
magnetic phase transitions discussed in Sec. 1111 B[ 

As can be inferred from Fig. [5l the resistivity measured along the a-axis is much smaller 
than p taken along the h- and c-axes. At 300 K, the ratio pb/ Pa ~ Pc/ Pa — 4. Moreover, the 
anomalies at Tn and Tt are hardly visible on the Pa{T) curve. This behavior is fully consistent 
with the established for CeRh3Si2 very large magnetocrystalline anisotropy, characterized 
by the magnetic moments confined to the he plane. 

Our several attempts to measure the resistivity of polycrystalline LaRh3Si2 failed because 
of multitude of microcracks present in the specimens. For this reason it was not possible 
to extract the magnetic contribution to the electrical resistivity of CeRh3Si2 employing a 
commonly used method of subtracting the resistivity of its non-magnetic counterpart. 



D. Specific heat 

Figure [6] displays the temperature dependencies of the specific heat of CeRh3Si2 and its 
non- magnetic isostructural counterpart LaRh3Si2. At low temperatures a distinct anomaly 
located at about 4.7 K confirms the antiferromagnetic ordering of the cerium moments in 
CeRh3Si2. The overall shape of C{T) obtained for the La-based compound is in turn typical 
for non-magnetic metals. In particular, below about 7 K it is easily describable by the 
formula: 

C-(T) = oT + ^||n (1) 

in which the first term is a conventional conduction-electron contribution to the specific heat 
with 7 being the Sommerfeld coefficient, and the second term is a low-temperature phonon 
contribution in a form of the Debye-T^ law with the characteristic Debye temperature Gn, 



the universal gas constant i?, and n being a number of atoms in the formula unit [34 1. 
Least-squares fitting of Eq. ([1]) to the experimental data (see the solid line in the inset to 
Fig. |6]) yielded the values of 7 = 3 mJ/(mol K^) and Gd = 428 K. At higher temperatures 
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FIG. 6: (Color online) (a) Specific heat of single-crystalline CeRh3Si2 (squares) and polycrystalline 
LaRh3Si2 (circles) as a function of temperature; the solid line is a fit of Eq. ([2]) to the experimental 
data, the arrow marks the Neel temperature. The upper inset shows C/T vs. measured for 
the La-based compound; the straight solid line is a fit of the Debye T^-law [Eq. ([I])]. The lower 
inset displays the specific heat of LaRh3Si2 divided by as a function of T; the arrow marks the 
position of a maximum (for its interpretation see the text). 



C {T) of LaRh3Si2 can be described by the equation: 



C{T) = 7T + no X 9i? ( 



;dx + nE X 3i? I — 



1)^ 



,2 ' 



(2) 



in which the second term is the full Debye expression for the phonon specific heat, and the 
third term describes the phonon specific heat within the Einstein model of lattice vibrations 



(see e.g. Refs. 



34 



J35| ) with Oe being the characteristic Einstein temperature. The coeffi- 
cients riD and n-^ are numbers of atoms in the formula unit (riD + ^E = n), vibration of which 
were assumed to follow the Debye or Einstein models, respectively. The solid line in Fig. [6] 
presents the temperature variation of the specific heat calculated employing Eq. ([2]) with the 
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FIG. 7: (Color online) (a) Logarithmic temperature variation of the magnetic contribution to 
the specific heat of CeRh3Si2 (left axis); the thin solid curve serves as a guide for the eye, the 
arrows mark the magnetic transition temperatures, and the thick solid line represents the Schottky 
contribution calculated for the CEF level scheme from Tab. IIVI The dashed curve is the magnetic 
entropy S of the compound (right axis); the drop lines mark the values of S{T^) and T{S = Rln 2). 
(b) AC/T at low temperatures; the solid curve is a fit of Eq. ^ to the experimental data in the 
magnetically ordered region, and the arrows mark the magnetic transition temperatures. 

parameters 7 = 3 mJ/(mol K^) and Od = 428 K (as determined from the low-temperature 
data), and Ge = 125 K, tzd = 5 and t^e = 1- The two latter values (riD and tie) were chosen 
arbitrarily, assuming that the rare-earth atoms (1/f.u.) in CeRh3Si2 and LaRh3Si2, which 
occupy the centers of the hexagonal tubes (see Fig. [1] and Sec. IIII A|) . are weakly coupled 
with their surroundings (5 at./f.u.) and exhibit some tendency to vibrate independently 
of the whole lattice. The characteristic temperature of those vibrations (i.e. the Einstein 
temperature) was estimated on the basis of simple analysis of C{T) measured for the La- 
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FIG. 8: (Color online) Specific heat of single-crystalline CeRh3Si2 measured at low temperatures 
in several different magnetic fields B applied along the 6-axis. The curves are shifted up for the 
sake of clarity. 

based compound: analytical properties of the Debye and Einstein functions allow to relate 
the position Tmax of the low-temperature maximum in the C{T)/T^ curve to the Einstein 
temperature as 9e ~ 5 x Tmax- In the case of LaRh3Si2 Tmax ~ 25 K, hence 9e ~ 125 K 
(see the lower inset to Fig. |6]). 

In order to estimate a 4/-derived contribution AC(T) to the total heat capacity of 
CeRh3Si2, the phonon specific heat of LaRh3Si2 was subtracted from C{T) of the Ce-based 
compound. As apparent from Fig. [7](a), the so-obtained AC{T) exhibits in the paramag- 
netic region a distinct Schottky-like anomaly due to CEF effect. The thick solid line in 
Fig. El^a) represents the Schottky specific heat calculated for the CEF model derived in 
Sec. IIIIEI Apparently, the model accounts reasonably well for the experimental data, ex- 
cept for high temperature region, where significant quantitative discrepancy is seen. The 
shortcomings likely arise because of the accumulation of experimental errors (see Fig. [6]) 
directly influencing reliability of the subtraction procedure. 

Below about 50 K, the magnetic contribution to the specific heat of CeRh3Si2 increases 
again, and at about 4.70(1) K AC{T) takes a form of distinct A-shaped anomaly, followed 
by another pronounced peak, that occurs at about 4.48(1) K. While the A-shaped feature 
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corresponds to the second-order antiferromagnetic phase transition evidenced at Tn in the 
magnetic susceptibihty data, the second feature can be ascribed to the spin-rearrangement 
at Tt. The characteristic spike-hke shape of that singularity is typical for first order phase 
transitions. In turn, the magnetic contribution to the specific heat in the range 5-50 K can 
be attributed to short-range interactions, which evolve finally into the long-range antifer- 
romagnetic ordering. The latter hypothesis is supported by the temperature dependence 
of the magnetic entropy, defined as S{T) = AC{T)/T dT. As can be inferred from 
Fig. [Tl^a) (right axis), the entropy at Tn is strongly reduced in comparison to the entropy 
of the ground state doublet (= i?ln2) and amounts only to about 0.44/2 In 2. The value of 
i?ln2 is achieved at temperature as high as ~ 30 K, suggesting that even well above the 
system is not yet completely disordered. 

In the ordered region, the specific heat of CeRh3Si2 is dominated by an antiferromagnetic 
magnon contribution. Assuming that the dispersion relation in the antiferromagnetic spin- 



wave spectrum has a relativistic form u = a/A^ -|- D'^k'^ (where u, D and k are the spin-wave 
frequency, stiffness, and wave number, respectively, and A is a spin-wave gap) AC (T) takes 



the form of a sum of the electronic and magnon contributions 36|] : 



AC(T) = 7T + AAV ^e-^/^ 



20 \T 32 VT 



(3) 



where A is a proportionality coefficient related to as A oc Least-squares fit- 

ting procedure [see the solid line in Fig. El^b)] yielded the values: 7 = 7 mJ/(mol K^), 
A = 3 mJ/(mol K^), and A = 16 K. The small value of the Sommerfeld coefficient, which 
can be also directly deduced from Fig. WO^) as AC{T)/T\^^q, is typical for simple metals. 

Figure E] shows the results of the specific heat measurements done in several different 
magnetic fields applied parallel to the 6-axis. As in the case of the magnetization mea- 
surements the latter experiments were also focused on the temperature region close to Tn 
and Tt. As seen, the two anomalies visible in zero magnetic field at 4.7 and 4.5 K split 
in higher fields into four separate features in C(T)/T, which independently shift to lower 
temperatures with rising B. However, all of the anomalies systematically broaden and in 
B ^ 0.7 T one can distinguish only two of them. Finally, in 5 = 1.2 T only one hump at 
about 3.7 K is visible, and in i? ^ 1.5 T the specific heat curve is featureless. The positions 
of the anomalies from Fig. |8] served as complementary set of data used for constructing the 
tentative phase diagram shown Sec. HVl 
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TABLE IV: CEF parameters, molecular field constant A, eigenenergies and eigenwavefunctions 
obtained from the fitting of the experimental magnetic susceptibility data and the Schottky con- 
tribution to the specific heat. 



Parameter 


Value [K] 








-B20 


1322 








B22 


-1265 








B40 


727 








B42 


1327 








B44 


-2345 








Bm 


-7478 








Bq2 


-2197 








Bq4 


3060 








Bm 


22239 








A [mol m"'^] 


501176 








SQX 


1.86% 








Eigenvalues [K] 


Eigenvectors 











- 0.801 15/2, 5/2) + 0.428 |5/2, 1/2) - 


■ 0.407 15/2, -3/2) 






655 


0.875 |5/2, 3/2) - 0.289 |5/2, -5/2) - 
... +0.239 |5/2, -1/2) - 


0.255 |7/2, 3/2)+ . . 
0.126 |7/2, -1/2) 






692 


- 0.783 |5/2, 1/2) - 0.464 |5/2, 5/2) - 


■ 0.334 |7/2, 1/2) + 0.223 |7/2, - 


7/2) 


3091 


- 0.765 |7/2, 7/2) - 0.511 17/2, -1/2) 


+ 0.313 |7/2, 3/2) - 


0.235 |7/2, 


-5/2) 


3478 


- 0.940 |7/2, 5/2) + 0.246 |7/2, 1/2) - 


■ 0.217 |7/2, -3/2) 






5467 


- 0.742 |7/2, -1/2) + 0.459 |7/2, 7/2) 


- 0.232 |7/2, 3/2) + 








... - 0.381 15/2, -1/2) - 


0.138 |7/2, -5/2) - 


0.133 |5/2, 


-5/2) 


5882 


- 0.854 |7/2, 3/2) - 0.383 |7/2, 7/2) - 


■ 0.233 |5/2, 3/2)+ . 








. . . + 0.202 |5/2, -5/2) + 0.170 |7/2, -5/2) 







E. CEF calculations 

The temperature dependence of the magnetic susceptibility (Sec. IIII Bp and the Schottky 
contribution to the heat capacity (Sec. IIIIDf) of CeRh3Si2 suggest predominantly locahzed 
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character of the /-electrons and the importance of definite sequence of the CEF levels. To 
capture the main features of the 4/^ state split in the CEF potential, we apply here the 
standard static description based on the renormalized Hamiltonian projected on the space 
of the pure /-electron states. The approach is known to be quite satisfactory for ins ulatin . 
systems, strongly supported by both qualitative and quantitative theoretical studies 



In the Wybourne notation 



37 



the one-electron CEF potential of the C2v point symmetry 



appropriate to the case of CeRh3Si2 reads : 



-ncEF — ^20^0 +-D22<-^2 " 



--D40t^0 " 



R 



R 



R r<(^)_LR 

(4) 



where Bkq are the CEF parameters, which are to be adjusted to the experimental data and 
are the normalized spherical harmonics. Besides, we employ also a simplified, angular- 



overlap-model (AOM) form of the CEF potentia' 



effective CEF interaction (see for instance Ref. 



38|): 



defined through matrix elements of the 



(m 



(5) 



rt fi 



In this equation, D^rmiO, Qrt, ^rt) is the matrix element of the irreducible representation D^"^^ 
of the rotation group, and Rrt, ©ri, ^rt are the angular coordinates of the nearest neighbor 
atom rt expressed in the coordination system set at the magnetic atom site. The index 
r distinguishes various chemically nonequivalent atoms and t runs over nearest neighbors 
r atoms. Thus, the whole CEF effect is described by three intrinsic parameters with 
fi = 0(0"), l(7r), 2(5) for each chemically nonequivalent nearest neighbor atom r. In the 
present case there are six silicon atoms and six rhodium atoms in the coordination sphere of 
the radius of 3.34 A. This coordination implies six non-equivalent AOM parameters, three for 
each nearest neighbor atom. If we neglect es and / or tie up with one third of (following 
the rules observed for insulating systems 37|), four-parameters or two-parameters versions 
of AOM can be applied. 

The Hamiltonian H includes also the spin-orbit interaction: 



so 



C4/I ■ s 



(6) 



with the spin-orbit constant (^4^- equal to 647.3 cm ^ (931.4 K) 40|, and the Zeeman term: 



M 



B 



B„ 



B. 



{Ix + 2s,) + (ly + 2Sy)^ + (/, + 2S,)^ 



B 



B 



B 



(7) 
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describing the system response to the external magnetic field B: 



H = Hqef + Hso + -^M 



(8) 



Following Refs. 



41 



and 



42 



the molar magnetization Mmoi,a and the magnetic susceptibility 



Xa are calculated from derivative of the free energy F of the system: 

F = 5^ exp (3 = l/{k^T) 

with respect to the magnetic field direction a, i.e.: 



M, 



mol,Q 



1 dF 
VdB 



where 



dEu.n 



dB 



(9) 



(10) 



(11) 



and Ey^a are the eigenenergies of the Hamiltonian given by Eq. |8l The molar susceptibility 
is defined as: 



M, 



mol,o 



B 



(12) 



The CEF parameters in Eq. H] have been fitted to restore the temperature dependences of 
the magnetic susceptibility measured along the a-, b- and c- axes and simultaneously account 
for the Schottky contribution to the heat capacity (see section lIllDI) . For this purpose the 
Levenberg-Marquardt method 43| was applied in the version implemented in the computer 



programs by Helmut Schilder and Heiko Lueken |41|J42|J44|J45| . In order to scan the entire 
parametric space, several initial CEF parameter sets have been probed. Among them one 
may distinguish three s pec ial cases: (i) zero-values, (ii) the CEF parameters reported for 



the compound CeIr3Si2 46|] and (iii) the parameters estimated using AOM with the values: 



-300, 300, 600 cm-i (-432, 432, 864 K), e 



Si 



-100, 0, 100 cm' 



-144, 



0, 144 K) and e^^ = -200, 200 cm'^ (-288, 288 K), e^'^ = ±^6^^, ef^ = -50, 0, 50 cm'^ 
(—72, 0, 72 K). Besides, several independent initial randomly generated CEF parameters 
sets (over 230) have been probed as well. 

The fitting procedure was same for each of the initial CEF parameters set. First, we 
fit the two second-order parameters, whereas the remaining ones were kept at their initial 
values. Then, the parameters of second and sixth rank were fixed at their values determined 
in the previous step, while the fourth-rank parameters were fit. Subsequently, all the second- 
rank and the fourth-rank parameters were varied simultaneously. In the following we call 
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this initial step the base phase, in which the general shapes of the experimental magnetic 
susceptibility curves have been reproduced fairly well, yet some details of the measured 
curves could not be accounted for. 

In the second, refining, phase we included the sixth-rank CEF parameters. Their influence 
on the thermal properties is expected to be of less significance, as they only determine the 
CEF admixture of the 7/2 states having energies above 2000 cm~^ (2878 K) to the low- 
energy 5/2 states predominantly governing the thermodynamic properties of the compound 
studied in moderate temperature ranges. This admixture modifies rather slightly both the 
CEF splitting of the 5/2 state and the matrix elements of the angular momentum operator. 
Nevertheless, the mixing may become more pronounced if the CEF effect is strong enough. 
Indeed, some preliminary simulations showed that in the case of CeRh3Si2 the admixture of 
the j = 7/2 states may reach as much as 10% for the lowest states. Hence, we decided to 



include the sixth-ran 



square error SQX [44 



i states in the refining phase, and as a result the relative root mean 
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48| has dropped from 2.92% to 1.86%. 



The fitting procedure of the sixth-order parameters were analogous to that of the fourth 
order. First, we varied only the sixth-order parameters, while keeping all the remaining 
parameters fixed on the values determined in the base phase. In the final step all the CEF 
parameters were varied simultaneously. In each step the molecular field constant A was also 
varied. 

Most of the probed initial sets of CEF parameters led to stable solutions of three different 
types. The corresponding fits yielded similar sequences of the CEF levels but distinctly 
different parameters sets, except for the two second-rank parameters, which appeared to 
have fairly similar values in each solution. The best fitting results are collected in Tab. IIVI 
As seen from this table, the CEF effect is generally large and leads to the total splitting 
of the 5/2 state of the order of 500 cm~^ (720 K) and the overall 4/^ spitting exceeding 
4000 cm~^ (5755 K). Despite the pronounced differences between the parameters in the 
three distinct sets (only the solution giving lowest SQX is shown) the sequence of levels 
and the model curves are almost same, what is reflected in negligible differences in the 
root-square errors and the energies. The scattering of some fourth- and especially sixth- 
rank CEF parameters in various solutions confirms rather limited reliability of the latter 
parameters. Consequently, as the final result we accepted the solution characterized by the 
lowest root-square error. The so-obtained model curves are shown in Fig. [2] together with 
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the experimental ones. Additionally, Fig. [Tl^a) presents the calculated Schottky contribution 
expected for the obtained CEF splitting scheme. 



F. X-ray photoemission spectroscopy 



Analysis of the Ce 3d XPS spectra is a proper method to investigate the character of 
the Ce 4/ states in Ce-intermetallic compounds due to the strong Coulomb interaction 
between photoemission holes in the 3d shell and electrons located near the Fermi level. 3d 
core ionization of Ce-based compounds usually results in a spectrum of_final states, whose 
configurations correspond to the /°, and p states (see Refs. 



49 



and 



50) 



Figure [9](a) shows the Ce 3d XPS spectrum obtained for CeRh3Si2. Two final-state 
contributions and are clearly observed, which exhibit spin-orbit splitting Aso — 18.6 
eV. In turn, the component is not evident in the spectrum, which implies stable valence 
3+ of the cerium atoms. According to the Gunnarsson and Schonhammer theory 49|], the 
hybridization energy Aj^, which describes the hybridization part of the Anderson impurity 



Hamiltonian 



5l| and is defined as TiV^pmax [where pmax is the maximum in the density of 



states (DOS) and V is the hybridization matrix element], can be estimated from the ratio 
r = I{p)/{I{f^)+I{p)) oi the I{p) and J(/^) intensities The particular contributions 
to the 3d XPS spectrum can be derived by its analysis in terms of the Doniach-Sunjic 
app roach 52| (the experimental details and accuracy of the method is also discussed in Ref. 
53l ). This way, for CeRh3Si2 one estimates a hybridization width Afs of ~56 meV. 

More evidence for the stable valence of Ce^"*" ions in CeRh3Si2 comes from the Ce 4d 
XPS spectrum displayed in Fig. Mjo). Clearly, there is no additional peak at hig her binding 
energies {E > 118 eV) that could be attributed to the Ce 3d^p state 
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54 



55| 



G. Electronic band structure calculations 



The band structure calculations were performed using the Wien2k code within the LDA, 
LDA-I-U and LSDA+OP schemes [the first two approaches with (LSDA and LSDA+U) and 
without spin polarization]. Non spin polarized densities of electronic states were used for 
calculations of photoemission spectra, which in experiment were measured at room tempera- 
ture, well above TAr=4.72 K. In the case of spin polarized calculations the antiferromagnetic 
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FIG. 9: (color online) (a) Ce 3d and (b) Ce Ad XPS spectra obtained for CeRh3Si2. The and 
components were separated on the basis of the Doniach-Sunjic theory. 

solution was not assumed in advance. The twelve atoms forming unit cell that was used in 
our calculations were treated as different types of atoms. The starting magnetic moments 
were chosen to have opposite sign. The system reached, iteration by iteration, a selfconsis- 
tent solution that was very close to the antiferromagnetic state. The resultant magnetization 
of both sublattices was almost balanced and was equal to about 2xl0~^ /i^ per f.u.. 

In the first step of our calculations, we employed the LSDA approach. The calculated 
densities of states (DOS) for CeRh3Si2 system are presented in FigHO) (a) the total DOS 
and (b) the contributions provided by 4f electrons located on cerium atoms forming two 
magnetic sublattices. The valence band can be divided into two parts below the Fermi 
level {^f)i here located at E=0. The first part located in the range of (-11.5;-7.5) eV is 
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formed mainly by electrons provided by silicon atoms Si(3s). The main part of the valence 
band between about -7.5 eV and Ep is created by Rh(4d) and Si(3p) electrons. The Ce(4f) 
electrons provide contribution nearest the Fermi level, and they have the main contribution 
to DOS at E=Ef: 5.43 states/(eV atom), 77.5% of the total value 7.01 states/(eV f.u.). 

Based on DOS at one can calculate the Sommerfeld coefficient 70, which in our case 
is equal to about 16.5 mJ/(mol K^). This value is much higher than the experimental one 
7exp=3 mJ/(mol K^), reported earlier in Sec. HID. This discrepancy is probably caused by 
inadequate description of the strongly correlated 4f electrons by the LSDA formalism. 

Hybridization changes initial configuration 6s^5d^4f^ to the final one 6s^'^^5d'^'^^4f^ '^^. To- 
tal number of electrons is different from those for pure cerium systems because of charge 
transfer from/to rhodium and silicon atoms as well as the charge accumulated in the inter- 
stitial region between atomic spheres. 

The values of site projected moments on cerium atoms are equal to +/-0.08 /i^ per Ce 
atom. The resultant value is not spectacular because of opposite alignment of spin and 
orbital moments. Until now the experimental values of local moments were not known. 

The shape of DOS obtained by LSDA+OP approach (Fig. (c) and (d)) is similar to 
that reported earlier but usually in this approach the magnetic moments increase, in our case 
to +/-0.5 fiB per Ce atom. Unfortunately the DOS at Ep is still too high: 10.65 states/(eV 
f.u.) with contribution of 4f electrons of about 8.52 states/(eV Ce atom). The Sommerfeld 
coefficient is equal to 25.1 mJ/(mol K^). 

On the two lowest panels of Fig. [10] we present results of LSDA+U approach with U 
equal to 6 eV. The additional U term to the LSDA approach shifts the 4f DOS (in LSDA 
located well below Ep) towards higher binding energies (now located around -2 eV). The 
unoccupied part of 4f DOS is moved to the energies about 4 eV above the Fermi level. In 
this way the contribution of Ce(4f) electrons to the DOS at EEp is practically negligible, 
and the Fermi level is located in a pseudo gap of the DOS. The value of DOS(Ei?)=1.01 
states/(eV f.u.) gives the Sommerfeld coefficient equal to 2.39 mJ/(mol K^). The latter 
value is slightly smaller than the experimental Sommerfeld coefficient as it should be for 
theory. The magnetic moments are equal to +/-0.09 fiB per Ce atom. 

Figure [TT] presents experimental and calculated X-ray photoelectron spectra based on non 



spin polarized DOS plots and calculated cross sections collected in Ref. 



28|. The best fitting 



to the experimental photoemission spectrum was obtained for the LDA calculations. The 
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FIG. 10: The total (per f.u.) and site projected (per Ce atom) densities of electronic states (DOS) 
for antiferromagnetic CeRli3Si2 calculated within the LSDA, LSDA+OP, and LSDA+U schemes. 



consistency between the calculated and experimental spectra is quite good. The spectrum 
is dominated by the Rh(4d) electrons. The Ce(4f) electrons provide contributions near the 
Fermi level. The peak located at about 8-9 eV below Ei? is formed mainly by the Si(3s) 
electrons (not shown here). 



IV. SUMMARY 



CeRh3Si2 was found to crystallize in the orthorhombic ErRh3Si2-type structure with 
honey-comb-like Rh-Si cages surrounding the Ce ions, which occupy just one inequivalent 
crystallographic position in the unit cell. The 4/ electrons are well localized and experience 
exceptionally strong (as for cerium intermetallics) crystalline-electric-field effect with sig- 
nificant admixture of two j'-multiplets (i.e. 5/2 and 7/2). In the whole temperature range 
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FIG. 11: The measured X-ray photoelectron spectrum for the valence band compared with the 
calculated ones within the LDA approach for CeR3hSi2 system. 

1 .4 1 1 1 1 1 1 1 




Temperature (K) 

FIG. 12: (Color online) Tentative B — T phase diagram for single-crystalline CeRh3Si2 based on 
anomalies observed in temperature dependencies of magnetic susceptibility (squares; Fig. [3]), field 
variations of magnetization (squares; Fig.|3|) and the temperature dependencies of the specific heat 
(circles; Fig. [8|) for B \\ 6-axis. Solid curves are hypothetical phase-boundary lines. Colors serves 
as a guide for the eye. 

covered CeRh3Si2 exhibits very large easy ah plane anisotropy in its magnetic and electrical 
transport behavior. At low temperatures the compound undergoes two subsequent anti- 
ferromagnetic phase transitions. The experiments performed in external magnetic fields B 
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applied parallel to the easy 6-axis revealed very complex B — T phase diagram of CeRh3Si2 
(Fig. [T2|) . which needs to be verified by neutron-scattering studies. The band structure cal- 
culations confirmed the antiferromagnetic ordering of the cerium magnetic moments. The 
calculated values of these moments are equal to 0.08, 0.09, and 0.5 /i^ per Ce atom for 
the LSDA, LSDA+U, and LSDA+OP approach, respectively. Only LSDA+U calculations 
yielded correct relation between theoretical [2.39 mJ/(mol K^)] and measured [3 mJ/(mol 
K^)] values of the Sommerfeld coefficients being less than one. The non-spin-polarized DOS 
plots reproduced well the experimental X-ray photoelectron spectrum of the valence band. 

It is worth noting that properties similar to those of CeRh3Si2 (i.e. "devilish" magne- 
tization, complex magnetic phase diagram, strong CEF effect) were reported for a number 
of magnetically ordered rare-earth intermetallics. Among them one can mention highly 
anisotropic CelrgSia TbNisGea ykd], TbNiaSis 1321 and CeSb |3o||3l|. All these com- 



pounds exhibit several successive field-induced metamagnetic phase transitions, which man- 
ifest themselves as a series of sudden jumps between well-defined plateaus of field-dependent 
magnetization. Detailed bulk direction-dependent physical properties measurements as well 
as neutron diffraction and x-ray resonant exchange scattering studies have revealed that 
in all of these compounds one deals with non-trivial field-induced spin rearrangements. In 
most cases the metamagnetic processes are relatively well understood and described phe- 
nomenologically, yet microscopic models were proposed only for a limited number of the 
systems (for survey see Ref. l57| and references therein). In most cases, the magnetic be- 
havior was interpreted in terms of the Ising model with next-nearest neighbour interactions 
(Axial Next-Nearest Neighbor Ising model, ANNNI), which (in an ideal case) can produce 
an infinite number of commensurate phases occurring with an infinite number of steps in 
field-dependent magnetization, called commonly a "devil's staircase" (for review see e.g . 



Ref. 



58l ). Other theoretical approaches were proposed for TmAgGe 59| and HoNi2B2C 60 |. 



The former compound has been found to exhibit some features of frustrated system with 
distorted antiferromagnetic Kagome-like layered structure. The sequence of field-induced 
metamagnetic transitions observed in the compound TmAgGe was successfully reproduced 
using a hamiltonian that describes only (i) relatively weak antiferromagnetic next-nearest- 
neighbor and ferromagnetic nearest-neighbor exchange adopted to the geometry of the sys- 
tem, and (ii) strong crystal electric field effect 6l|. In turn, HoNi2B2C is a magnetically 
ordered superconductor, which shows at low-temperatures a transition from commensurate 
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antiferromagnetic to incommensurate c-axis complex spiral state with some a-axis incom- 



mensurate modulation (see Refs. 



62-64 



line 



and references therein). Here again the crysta 
electric field was taken into account to explain the complex magnetic phase diagram 
yet the generalized ANNNl model (a kind of so-called clock model) was also successfully 
applied [66 1. 

The microscopic origin of the complex magnetic behavior of CeRh3Si2 cannot be defined 
at this stage of our study. The most relevant factor governing its strongly anisotropic 
properties in the ordered and paramagnetic states is clearly unusually strong crystalline 
electric field potential that yield the total splitting of the 5/2 state of about 720 K and the 
overall 4/' splitting as large as about 5760 K. The estimated magnitude of the CEF effect in 



ics. The other cerium 
and CeIr3Si2 
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CeRh3Si2 is perhaps the largest ever reported for Ce-based intermetal 
compounds known for their giant CEF interactions are CeRh3B2 67 
where the total CEF splitting is about 2000 K and 1600 K, respectively. Remarkably, the 
latter compound is isostructural to CeRh3Si2, while the crystal structure of CeRh3B2 is 
closely related to that of CeRh3Si2. Further studies, including inelastic neutron scattering 
experiments and ANNNI calculations, are needed to shed more light on the unique features 
(i.e. unusually strong CEF and complex magnetic behavior) of the latter compound. 
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